Self-isospectrality, special supersymmetry, and their effect on the band structure 
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We study a planar model of a non-relativistic electron in periodic magnetic and electric fields that 
produce a ID crystal for two spin components separated by a half-period spacing. We fit the fields 
to create a self-isospectral pair of finite-gap associated Lame equations shifted for a half-period, and 
show that the system obtained is characterized by a new type of supersymmetry. It is a special 
nonlinear supersymmetry generated by three commuting integrals of motion, related to the parity- 
odd operator of the associated Lax pair, that coherently reflects the band structure and all its 
peculiarities. In the infinite period limit it provides an unusual picture of supersymmetry breaking. 
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Introduction - Supersymmetry as a fundamental sym- 
metry of Nature still waits for experimental confirmation, 
but as a kind of symmetry between bosonic and fermionic 
states, it already turned out to be fruitful in diverse areas, 
including nuclear [J, atomic, solid-state, and statisti- 
cal physics Q. Supersymmetric quantum mechanics 
was introduced under investigation of the problem of su- 
persymmetry breaking in field theory Q . In the simplest 
case a system is characterized there by a 2 x 2 diagonal 
matrix Hamiltonian, H = diag{H'^ , H^), and by two 
anti-diagonal matrix integrals of motion (supercharges) 
Qi and Q2 — ia^Qi. Supercharges are first order differ- 
ential operators generating an iV = 2 superalgebra 



(1) 



Such a system has an additional integral of motion F = 
(73, — 1, which classifies the states with F = +1 and 
— 1, by convention, as bosonic and fermionic states. Since 
[r, H] — and {F, Qa} — 0, the Hamiltonian and super- 
charges arc identified as bosonic and fermionic genera- 
tors. The cases of unbroken and broken supersymme- 
try are distinguished by the Witten index A defined as 
the difference between the total numbers of bosonic and 
fermionic states. In a non-periodic one-dimensional sys- 
tem, unbroken supersymmetry is characterized by one 
singlet ground state of zero energy and A 7^ 0; for bro- 
ken supersymmetry there is no zero energy singlet state 
and A = 0. However, as it was observed in [1, su- 
persymmetric periodic systems may support two zero en- 
ergy ground states, and then A = in the unbroken 
case. Some systems investigated there possess a specific 
property of self-isospectrality, meaning that correspond- 
ing superpartner potentials V'^ and V~ are given by the 
same periodic function shifted for a half of a period 2L, 
V~^{x-\-L) = V~{x). They belong to a class of finite-gap 
periodic systems, which play an important role, in par- 
ticular, in condensed matter physics and in the theory 
of nonlinear integrable systems Q . 

In a non/meargeneralization of supersymmetric quan- 
tum mechanics supercharges Qa are higher (n > 1) 



order differential operators generating a nonlinear super- 
algebra {Qa,Qb} = 2SabPn{H), With Pn{H) a polyno- 
mial of order n. The number of singlet states can take 
there any value from to n, and, as in periodic models 
with linear supersymmetry, the Witten index does not 
characterize supersymmetry breaking 0, . This indi- 
cates that in periodic finite-gap systems nonlinear super- 
symmetry may play an important role. 

To investigate the question of the presence and nature 
of nonlinear supersymmetry in periodic finite-gap sys- 
tems, in this Letter we study a planar model described 
by the Pauli Hamiltonian for a non-relativistic electron 
in periodic electric and magnetic fields. The model be- 
longs to a broad class of periodic systems investigated by 
Novikov et al. [11]. It is well known that in the absence 
of an electric field the model, which includes the Landau 
problem as a particular case, is characterized by a super- 
symmetry with the usual linear superalgebraic structure 
(|T]) Q . We choose periodic magnetic and electric fields in 
such a form that the spin-up and -down components of 
the electron wave function feel the same one-dimensional 
effective periodic potential but with a shift of half of the 
period. As a result, the effective potential of superex- 
tended system satisfies a property of self-isospectrality. 
Vector and scalar potentials are fitted to produce the as- 
sociated Lame equation with two integer parameters m 
and I, which belongs to a broad class of finite-gap systems 



with a smooth potential, see Eq. ([3|) below (8|, 12, [T^ . 

We find here a special nonlinear supersymmetry of the 
previously unknown structure, in which all the peculiar- 
ities of the band structure of the system are imprinted. 
In the infinite-period limit our system provides an un- 
usual picture of supersymmetry breaking rooted in its 
nonlinearity. 

Model - Consider a non-relativistic electron confined 
to a plane and moving in the presence of an electric field, 
given by a scalar potential (j){x,y), and a perpendicu- 
lar magnetic field Bz{x,y). It is described by the Pauli 
Hamiltonian 



He^{p.:, + A^) +{py + Ay) +aj.B,^. 



(2) 
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where the units are h = c = 2m = — e = 1. Let us restrict 
Bz and (p by the condition that they depend only on x. 
We choose = 0, Ay ~ w{x), then = Present 
the wave function in the form '^{x,y) = e^'^yip{x), where 
K, —oo < K < oo, is the eigenvalue of Py. Taking 
w{x) — a^ln(dnx) and (f>{x) = [3uP'{x) + ^w{x) + 5, 
with appropriate choice of constant parameters a, f3, 7 
and S, we reduce ([2]) to a quantum periodic system given 
by the diagonal matrix Hamiltonian H with up (+) and 
down (— ) components of the form H^^ j = — 3^2 
RereV+,{x) = V-,{x + L), 



'J 2 



-Cmdn X - d 



J 2 
dn X 



(3) 



Cm — m{m + 1), Ci = l{l + 1), m and I are integers 
such that + Cf 7^ 0, c is a constant; dux = dn (x, k) 
is the Jacobi even elliptic function, satisfying a relation 
dn (x + K) — fc'/dnx, with modulus < fc < 1 and 
real and imaginary periods 2K and AiK' , K{k) is the 
elliptic complete integral of the first kind, K' = K{k'), 
and k' is a complementary modulus, k'^ = 1 — k^. The 
Hamiltonian H obtained in this way then describes a pair 
of two parity- even associated Lame systems, shifted one 
with respect to the other for the half of the real period 
2L of the potential, that is equal to 2K for Cm 7^ Ci , and 
K for Cm = Ci. The cases with CmCi = correspond 
to the Lame system Q. By the Landen transformation 
14| , the case Cm = Ci is reduced to the case of the Lame 
system with same value of parameter Cm but with Ci — 
0, and with the same imaginary period but the real period 
divided in two. Taking into account that C-m-i = Cm, 
without loss of generality one can assume that m > I > 0. 
Isospectral subsystems ^ and H~ ^ belong to a class 
of finite-gap systems with even potential and number of 
energy gaps in the spectrum equal to m [1, . 

Hidden bosonized supersymmetry - Consider an n-gap 
periodic system with even Hamiltonian H = —-^^ + 
V{x), V{x + 2L) = V{x) = V{~x). Its spectrum a{H) 
is characterized by the band structure cr[H) = [Ea,Ei\ U 

. . . U [E2n-2,E2n-l\ U [E2n,0o), Eq < Ei < . . . < E2n, 

which consists of n valence bands and a conduction band, 
separated by energy gaps corresponding to n prohibited 
bands. The 2n + 1 singlet band-edge states of definite 
parity and energies Ej, j = 0, . . . , 2n, are given by pe- 
riodic or antiperiodic states of periods 2L or 4L. The 
states in the interior of permitted bands are described by 
Bloch-Floquet quasi-periodic functions, and every inter- 
nal energy level is doubly degenerate. The energy dou- 
blets are distinguished by the reflection (parity) operator 
R, Rip{x) — ip{~x), that is a nonlocal integral of motion. 
On the other hand, double degeneration of energy levels 
is a characteristic feature of a quantum mechanical N = 2 
supersymmetric system. The presence of 2n-fl singlets is 
an indication of the higher order (> 3) nature of the cor- 
responding hidden supersymmetry. Any finite-gap sys- 



tem is characterized by a nontrivial integral of motion 
A2n+i that is a self-conjugate differential operator of the 
formy42„+i ^i ^J„+i + a2n-iix) ^J„_i +...ao{x). The 
{A2n+i, H) is the Lax pair of the n-th order Korteweg-de 
Vries (KdV) equation, and the condition [A2n+i,H] = 
defines the stationary KdV hierarchy. This pair of com- 
muting operators satisfies identically the relation Q 

2n 

^2„+l=^2„+lW, P2n+l{H) = Y{{H-E,), (4) 



and 2n + I singlet states are the common eigen- 
states of H and A2n+i of the eigenvalues Ej and 0. 
The square of the self-conjugate operator A2n+i is pos- 
itive semi-definite, and ([4]) implies the described band 
structure of the spectrum. Integral ^2n+i is parity-odd, 
{i?, A2n+i} = 0. Taking into account that [R, H] = and 
i?^ = 1, we see that here the reflection operator plays the 
same role as the operator (73 for N = 2 superextended 
matrix system, and the operator Z — A2n+i can be iden- 
tified as the supercharge. Define a nonlocal odd operator 
Z = iRZ. Odd supercharges Qi = Z , Q2 = Z generate 
the N = 2 nonlinear supersymmetry of order 2?! -f- 1 : 



{Qa,Qb} = 2SabP2n+l{H). 



(5) 



Since this structure appears in the one-dimensional sys- 
tem without matrix (spin) degrees of freedom, the de- 
scribed nonlinear supersymmetry of order 2n + 1 for 
any n-gap periodic system with parity-even Hamilto- 
nian is identified as a hidden bosonized supersymmetry 

Sin US- 

Supersymmetry and band structure - Let us return 
to our finite-gap self-isospectral system. A vector space 
spanned by singlet band-edge states of the subsystem 
I or I is divided into two vector subspaces formed 
by 2L-periodic and 2L-antiperiodic (i.e., 4L-periodic) 
states. The singlet state Vl/o with the lowest energy Eq 
is 2L-periodic, and the singlet state of the other edge 
of the first valence band, ^^i, is antiperiodic. The two 
edge-states ^'2^-1 and ^2j, j = 1, ■ • . , m, separated by 
an energy gap, have the same period @, [I4I. Therefore, 
the space of periodic singlet states has odd dimension, 
and the space of antiperiodic states has nonzero even di- 
mension. On these two subspaces of singlet states, two 
irreducible non- unitary representations of the sl(2, M) al- 
gebra are realized. Namely, according to the 
space of 2m + 1 singlet states of the associated Lame 
system with m > I can be treated as a direct sum of 
two s^(2, R)-representations of dimensions m — I (spin 
ji — ^{m — I - 1)) and m + I + 1 (spin j2 = ^{m + I)). 
The period of the states of these subspaces is dictated 
by the parity of m — I. When m — Ms odd, spin-ji 
(spin-j2) representation is realized on 2L-periodic (4L- 
periodic) states, for even m — I the periodicity of spin-ji 
and spin-j2 subspaces interchanges. Making use of the 
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two corresponding algebraization schemes [l8|, we find 
two commuting anti-diagonal self-conjugate integrals of 
motion, X and F, where 



m - / = 4 •[ 

1 



Y = 1" 











(6) 
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cn 



cn X d 



i\x dx 



m+l + 1 



dn' X \^+ ( \ 



Y~i{x + K), ey = 1 (0) for m + I even (odd), and X 
has a form similar to ^ with -'^'„;(a;) = 
X+ i{x) = X^^ iix + K), ex = I ~ ey. The order of the 
differential operator X, \X\ = ni — l, is less than the order 
of Y , \Y\ = \X\ + 21 + 1 = m + I + 1. X and Y have op- 
posite parities (-1)™-' (X), and (-1)™-'+! (Y). When 
m — I is odd, 2L-periodic (4L-periodic) singlet states of 
each subsystem are zero modes of the integral X (Y). For 
even m — I the role of the operators X and Y as annihila- 
tors of periodic and antiperiodic edge states interchanges. 
The anticommutator {X, Y} — 2Z produces a diagonal 
integral Z = dza5(Z+ ;,Z„ ^), \Z\ - \X\ + \Y\ = 2m+l, 
with down, 



T7l , / V 



m,l 



and up, i{x) = Z^i{x + K), components, which 
are the parity-odd integrals of motion annihilating all 
the singlet states of the m-gap subsystems -ff ~ ^ and 
^mi described above. Hence, Z"^ = Pz{H), where 
Pz (H) is a spectral polynomial ^ of order 2m + 1 with 
H = diag {H^ H~ From the exphcit form of the in- 
tegrals X, Y and Z one finds that [Z,X] = [Z,Y] = 0, 
and X^ = PxiH), Y^ = Py{H). The polynomials 
Px{Pt) and Py{H) factorize the spectral polynomial 
Pz{H), Pz{H) = Px{H)Py{H), and include those fac- 
tors {H — Ej) for which Ej^s are the eigenvalues of band- 
edge states of corresponding periodicity [l^l . 

Periodic (antiperiodic) states have an even (odd) num- 
ber of nodes in the period interval. The maximal num- 
ber of nodes that can have the band-edge states an- 
nihilated by X^ I and Y^ ^ is not more than the or- 
der of these differential operators. When m — I \s 
even, X annihilates m — I anti-periodic band-edge states 
with 1, 1, 3, 3, . . . , m — / — l,m — I — 1 nodes. The 
m + / + 1 periodic edge states annihilated by Y have 
0, 2, 2 . . . , (m-|-Z), (m-|-Z) nodes. We find that when ra — l 
is odd, X annihilates m — / periodic band-edge states with 
0, 2, 2, . . . , m — / — 1, m — / — 1 nodes ii m — I > 1, and 
one nodeless state '^q if m — Z = 1. In the last case 
X is the usual first order supercharge [13]. The oper- 
ator Y annihilates m + l + 1 anti-periodic states with 
1,1, ... ,171 + l,m -\- 1 nodes. 

The picture can be summarized as follows. The band- 
edge state ^0 is a zero mode of the parity-odd super- 
charge, i.e. of Y {X) when m — Z is even (odd). The 
band-edge states of the same prohibited band 'attract' 
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FIG. 1: (color online). Scheme of band structure for self- 
isospectral systems with m = 4. Triangles (dots) indicate 
band-edge states annihilated by X (Y), the digits below mean 
their node numbers. The states with even (odd) number of 
nodes are periodic (anti-periodic). 



each other; they appear as zero modes of the same su- 
percharge. When the number of permitted bands m+1 is 
fixed, and m — I increases in steps of 2, there appear two 
new band-edge states annihilated by X, with increasing 
energies. Every such pair is separated by a pair of zero 
modes of Y . The highest 2{l -\- 1) singlet states are zero 
modes of Y . These properties are illustrated on Fig. [TJ 

Integral Z refiects the degeneration of the states of each 
subsystem, while X and Y reveal the self-isospectrality 
of the composed system. As a result, it is characterized 
by the 4-fold degeneration of quasi-periodic states and 
the double degeneration of the band-edge states. 

Nonlinear superalgebra - Besides nontrivial integrals 
X, Y and Z, our system is characterized also by mutu- 
ally commuting integrals Ti = a3,T2 = R and Fa = a^R. 
Any of them can be chosen as the operator F that clas- 
sifies all the integrals into bosonic and fermionic opera- 
tors. Appropriate linear combinations of physical states 
for which F = +1 and —1 are identified as bosonic and 
fermionic states. Let us choose F = 0-3. Eight integrals 
X, TiX, Y and TiY, i — 1,2,3, anticommuting with F 
are identified as fermionic operators. They anticommute 
between themselves for certain linear combinations of the 
8 bosonic operators Z, TiZ, Ti and H with coefficients 
that are some polynomials in H. Linear combinations 
of the bosonic operators J^i^^ = —^Ra^ZH^, ^jj^'' = 
icr3^n± and Jg^^^ = -\Rn±, where n± = i (1 ± era), 
have the only nontrivial commutators 



(8) 



Here a,h,c = 1,2,3, /9i,2 — 1, pz — Pz{H). This is a 
nonlinear deformation of su{2) © su{2) © u{l) © u{l), 
where the two last terms correspond to F = CT3 and 
H. The nonlinear algebra ([8]) is reminiscent of nonlinear 
symmetry algebra generated by the angular momentum 
and Laplace-Runge-Lenz vector operators in the quan- 
tum Kepler problem 21|. The complete superalgebra 



is identified as a nonlinear deformation of the su(2|2) 
superunitary symmetry, in which H plays a role of the 
multiplicative central charge (2^ . 

Infinite period limit - In the self-isospectral system 
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FIG. 2: (color online). Qualitative picture of supersymmetry 
breaking in a self-isospectral system with m = 3, / = f in 
the inhnite-period limit. The form of potentials and bands 
shown on the left corresponds to the modulus k close to 1. 
Two lower horizontal dashed lines on the right show energy 
levels of singlet bound states, the upper separated horizontal 
continuous line corresponds to a doublet of bound states, the 
line at the bottom of continuous spectrum indicates a doublet 
of the lowest states of the scattering sector. 



considered here, band-edge states form energy doublets, 
quasi-periodic states are organized in quadruplets. In 
the infinite period limit, corresponding to fc — > 1, fc' — > 
0, K oo, dn{x,k) — > ^^^^ ^ , the system transforms 
into a pair of reflectionless Poschl- Teller systems given by 
potentials V^^ = -C=^cosh"^a; c, C+ = C;, C" = C™. 
In this limit the spectral polynomial degenerates [l^. 
The quasiperiodic states of the conduction band reduce 
to the scattering states, supersymmetric doublet of its 
band-edge reduces to a doublet of the lowest states of the 
scattering sector. In a shifted subsystem ^ , m—l lower 
valence bands disappear, while the rest of them in both 
subsystems shrink to the bound states. The resulting 
system is characterized by m — Z singlet and I doublet 
bound states, and by a doublet of the lowest states of 
the scattering sector. The rest of the scattering states is 
organized in energy quadruplets. This unusual picture of 
supersymmetry breaking, illustrated on Fig. O is related 
to the nonlinear nature of self-isospectral supersymmetry. 

Conclusion - In the model investigated here the non- 
linear self-isospectral supersymmetry originates from: i) 
the separability of the singlet band-edge states of both 
subsystems into two non-empty subspaces of periodic and 
anti-periodic states, and ii) the related factorization of 
the higher order Lax pair operator of the associated sta- 
tionary KdV hierarchy. The mutually commuting in- 
tegrals X and Y are the annihilators of the band-edge 
states of definite periodicity. They factorize the integral 
Z that annihilates all the band-edge states. The unusual 
nonlinear supersymmetry generated by these nontrivial 
integrals together with integral and parity operator R, 
reveals the band structure of the system and all its pe- 
culiarities in the same way as a the nonlinear symmetry 
associated with the Laplace-Runge-Lenz vector reflects 
specific properties of the hydrogen atom spectrum 21|. 
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